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We study the point of transition between complete and incomplete financial models 
thanks to Dirichlet Forms methods. We apply recent techniques, developped by Bouleau, to 
hedging procedures in order to perturbate parameters and stochastic processes, in the case 
of a volatility parameter fixed but uncertain for traders; we call this model Perturbed Black 
CN ■ Scholes (PBS) Model. We show that this model can reproduce at the same time a smile 

effect and a bid-ask spread; we exhibit the volatility function associated to the local- volatility 
fV| model equivalent to PBS model when vanilla options are concerned. Lastly, we present a 

f-H ■ connection between Error Theory using Dirichlet Forms and Utility Function Theory. 

£H , Key Words: Smiles, dynamic hedging, local volatility, stochastic volatility, error, Dirich- 

let form, carre du champs operator, bias, utility functions. 

1 Introduction 

> : 

r**- . In this article, we study the impact of uncertainty of volatility on the price of vanilla options. 

A classical approach consists in perturbating volatility by a small random variable and in 
O ' considering the associated truncated expansion, but this approach encounters difficulties in dealing 
\0 ' with infinite dimensions. An alternative way, based on Dirichlet Forms (for reference see Albeverio 
P, Bouleau et al. [3] and Fukushima et al. [12]), has been suggested by Bouleau [I]; this 
O ■ idea yields the same representation of small perturbations as that of the classical perturbation 
j> \ approach. 

In classical theory of financial mathematics, we assume that all market securities have a definite 
price. Indeed, the hypothesis of completeness of the market (see Lamberton et al. pj)]) forces a 
single price for a contingent claim. If we take into account an uncertainty on a parameter, we find 
that the price of the contingent claim is not unique but we have many possible prices, therefore 
we can reproduce the bid-ask spread by means of a utility function related to the uncertainty on 
prices. 

If the uncertainty on parameter is small, we may neglect orders higher than the second, so we 
chose to work with Gaussian distributions. In Error Theory using Dirichlet Forms, we associate 
the variance to the "carre du champ" operator T and the shift to the generator of semigroup A. 

Historical Black Scholes model for asset pricing assumes that the diffusion process for asset 
price is log-normal with a constant volatility; however, many works on empirical market data 
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present a skewed structure of market implied volatilities with respect to the strike; this effect is 
called smile of volatility or volatility skew (see Renault et al. [18] and Rubinstein [19]). 

Implied volatility is convex as a function of the strike and generally exhibits a slope with 
respect to the strike at forward money (see Perignon et al. [T7]); to take this into account, we 
propose a new model based on BS model, characterized by an uncertain volatility parameter, 
called Perturbed Black-Scholes model. This is a stochastic volatility model with closed forms for 
option pricing; we study some constraints to force a smile on implied volatility and define the 
local volatility model, by means of its volatility function, equivalent to the PBS model for vanilla 
options. 

Summarizing, we propose a new financial model for securities pricing based on the Black 
Scholes model with a random variable as volatility, we use a perturbative approach to preserve 
closed forms for options prices and greeks; this model permits to reproduce a smile on implied 
volatility and generate automatically a bid-ask spread. 

The paper is organized as follows: 

In section 2 and 3, we present two ways to perturb the volatility, the classical approach and 
the one using Dirichlet Forms. In section 4, we present the PBS model and study the effect of 
uncertainty on volatility for an underlying following Black Scholes model without drift. In section 
5, we investigate the relations with the literature, while in section 6, we present an interpretation 
of the relative index defined in section 4 by means of utility functions theory. Finally section 7 
resumes and concludes. 

2 Classical Approach 

We consider a function F, which denotes the payoff or the hedging value of an option and depends 
on a parameter a (the volatility for instance), then we consider a perturbation of the parameter 
by means of a "normal" distribution and analyze the impact of this perturbation on F; we suppose 
that F € C 2 with respect to the parameter a. 

This approach may be performed with elementary limit calculations in a finite dimensional 
framework: 

One-dimensional Case 

Assume that a £ R. We model the perturbation on the parameter a by means of the following 
transformation: 

cro -> <T + yjeyg + ea 

where g is a random variable, following the centered reduced normal law and independent of all 
the other random variables present in this financial problem, like the Brownian Motion; a is a 
constant and e is a very small parameter. In other words we replace the parameter with a random 
variable with mean <7o + ea and variance ej. We can interpret the term 7 as the normalized 
variance of the variable a, subject to the perturbation; a is the bias of the parameter induced by 
the perturbation. 

We evaluate the bias and the variance of the pertubation of F, i.e. F(a + ^JT^g + ea) — F(ao). 
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Finally we can remark: 
Remark 2.1 T7ie &ias and variance have the two following chain rules: 

1. the bias of the function F, induced by the parameter uncertainty, depends both on the bias 
and the variance of the parameter. Indeed, this bias presents two terms. The first term, 
given by the bias of the perturbed parameter, is proportional to the first derivative. The 
second term is related to the convexity of the function F and proportional to the variance of 
the parameter; it has a purely probabilistic origin (see Bouleau ^ and JE/). 

2. the variance of the function F is proportional to the variance of the perturbed parameter and 
to the first derivative of the function. 



The figure [H resumes the two impacts. 




Figure 1: Impact of uncertainty on a parameter through a non-linear function. 
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Multi-dimensional Case 

We can extend the previous results to the case of a ^-dimensional parameter; assume a"o G M. d 
and F : M. d — > R. In this case the part of the parameter is expressed as 



where 



G ~ jV(0, I d ) 

7 G M(d) symmetric and positive definite 
72 is a square root of 7 
a G R d 

We suppose again that e G R is a small parameter, so that we can approximate the bias and 
the variance 
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where 7^ ; , indicate the element (k I) of matrix 72 , A, indicate the i t/l component of vector A and 
likewise Gj is the i th row of random matrix G. 

Finally we find a similar result in one- dimensional case: 

Remark 2.2 The bias and variance have the two following chain rules: 

1. the bias of the function F induced by the uncertainty of the parameter depends both on the 
bias and the variance of the parameter, as a matter of fact this bias presents two terms: a 
first term given by the scalar product of the perturbed parameter bias and the gradient of 
function F; the second term is related to the convexity of the function F and proportional to 
the covariance of the parameter a; 
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2. the variance of the function F induced by the uncertainty of the parameter is proportional 
to the norm of gradient of F weighted by the covariance matrix of the perturbed parameter 
a. 

Several ways are conceivable in order to extend this approach to infinite dimentional parame- 
ters. This is compulsory since we are generally concerned with perturbation of a stochastic process. 
One way is based on Dirichlet Forms which possess several convenient featurs for computing the 
propagation of perturbations. 

3 Dirichlet Forms Approach 

Then we present the Dirichlet Forms approach, we start by recalling the essential ingredients of 
this method (see Bouleau [1]). 
We define an error structure: 

Definition 3.1 (Error structure) An error structure is a term 

(p, p, d, 

where 

Q, J 7 , pj is a probability space; 

D is a dense sub-vector space of L 2 J 7 , P 



• T is a positive symmetric bilinear application from D x D into L 1 [Q, J 7 , PJ satisfying the 

functional calculus of class C 1 fl Lip, i.e. if F and G are of class C 1 and Lipschitzian, u and 
v G D, we have F(u) and G(v) G D and 

T [F(u), G(v)} = F'(u)G'(v)T[u, v] P a.s.; 

• the bilinear form £[u, v] = ^E\T[u, v}} is closed; 

• The constant function 1 belongs to D, i.e. the error structure is Markovian. 

In mathematical literature the form £ is known as a "local Dirichlet form" that possesses a 
"carre du champ" operator Y. 

We recall the definition of sharp operator associated with V 

Definition 3.2 (Sharp operator) Let {vt, T . P, D, rj an error structure and (q, J 7 , Pj a 

copy of the probability space (p, J 7 , Fj . Under the Mokobodzki hypothesis that the space D is 
separable, there exists an operator sharp ( )* with these three properties: 

• VueD, m # 6 L 2 (P x P); 

• Vm G D, T[u] = E \(u*) 

• Vu G W and F G C 1 n Lip, (F{ Ul , ... , u n ))* = J2ti (g ° u ) u t ■ 



From a numerical point of view, the sharp operator is an useful tool to compute T because the 
sharp is linear whereas the carre du champ is bilinear. We can also associate to the error structure 

\Q, J 7 , P, D, a unique strongly continuous contraction semi-group (Pt)t>o y i a Hille Yosida 
theorem (see Albeverio p] pages: 9-11 and 20-26). This semigroup has a generator (A, DA); it 
is a self-adjoint operator that satisfies, for F £ C 2 , u £ DA and T[u] £ L 2 (P): 

A [F(u)] = F'(u)A[u] + -F"(u)T[u] P a.s. 

The functional calculus extends the ideas of classical Gauss error theory, the idea is to consider 
the perturbation as an error. In analogy with the classical approach of error theory we associate 
the carre du champ operator T to the normalized^ variance of the error, the sharp operator 
becomes a linear version of the standard deviation of the error. Similarly, the generator describes 
the error biases after normalization (for more details we refers to the book of Bouleau |lj chapters 
III and V or to [5]). 

3.1 Perturbation on a parameter 

Now we go back to the study of the impact of a perturbation. To do so we implement an error 
structure on the parameter. Therefore we consider the same function F £ C 2 depending on a 
parameter, the volatility for instance. We model the perturbation on the parameter with the 
following transformation: 

o"o — > o = O"o + X. 

Here X is a random variable that represents the error on the parameter a. We consider an 
error structure associated to X with some hypotheses to allow the computations. 
Assumptions: 

1. X £ D, T[X] and A[X] are known; 

2. functions x i— > T[X](x) and x t— > A[X](x) belong to L 1 H L 2 are continuous at zero and not 
vanishing; 

3. the error structure possesses a gradient operator, allowing to define the sharp operator which 
is a special case of gradient (see Bouleau et al. [3] chapter V § 5.2 or Bouleau [1] chapter V 
§2). 

Denoting A* the sharp of the parameter a, we have 

F(a)* = F\a)X* 
T[F(a)} = [F'(a)] 2 T[A] 

A[F(a)] = F>{a)A[X] + l -F"{a)T[X]. 

Following the idea of truncated development, we are interested in evaluating the bias and the 
variance at the value a = (To; i.e. at X = 0: 

T[F(a)](a ) = [F'(a )} 2 T[X]\ X=0 
A[F(a)](a ) = F'(a )A[X]\ x=0 + ^F"(a )T[X]\ x=0 . 

1 the variance divided by the small parameter e. 
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We can remark that we have found the same chain rules for bias and variance as in the classical 
approach; however it is worth notice that we had not to specify the dimension of the space where 
X is defined. This dimension can be finite or infinite, as in thecase of the Wiener space, where we 
are able to define coherent error structures (See Bouleau et al. [3] and [3]). Hence this framework 
is an extension of the classical approach. 

Now we present a classical case where we can perform explicit computation. 
Example 3.1 (Diffusion) 

Proposition 3.1 Let X t be a diffusion process of the form 

dX t = a{X t )dB t + b{X t )dt 

where a(.) and b(.) satisfy classical conditions and B t is a brownian motion. We use this SDE to 
perturb a parameter X , the time t play the role of the small parameter e and we are interested in 
studying the impact of this perturbation on a smooth function F depending on the parameter X , 
we have the following bias and variance: 



1 
2 

\X=X = i r W)l a ^o)- 

The bias is given by the generator of the diffusion and the variance by the martingale term. 



A[F(X)]\ x=Xq = F>(X )b(X ) + -F"(X o y(X 0) 
nF(X)]\x=x n = [F'(X )} 2 a 2 (X 



We use a short lemma: 
Lemma 3.1 Let M t a stochastic integral: 



-t 

Mt = I K«LJB« 



t I s tJ ' ± - > s 



where K t and L t are continuous squared integrable processes and N t the following stochastic scaled 
brownian motion: 

N t = K L B t . 

Then N t is an approximation of M t in L 2 when t goes to zero. 
Proof of lemma: Indeed 



jo) 2 ds 



[(M t - N t ) 2 ] = - f (K S L S - K L C 
t t Jo 

< - [ {[K s ] 2 [L s - L ] 2 + [K s - K } 2 [Lo] 2 } ds - □ 
t Jo 

Proof of proposition 3.1: By Ito formula 



F(X t ) = F(Xo) + / F'(X s )a(X s )dB s - 
J o 

= F(X ) + [ A [F(X S )} ds + Mi 
Jo 



F'(X s )b(X s ) + l -F"(X s )a 2 (X s ) 



ds 



where M t is a martingale. 

Thanks to the previous lemma, we approximate M t with N t , given by 



N t = F'(X )a(X )B t 

Finally we can write the approximate law of F(X) generated by the perturbation on the 
parameter X 

F(X t ) « F(X ) + tA[F(X )\ + VtF'(X )a(X )B 1 
where F(X ), A[F(X )] and F'(X Q )a(X Q ) are .Fo-measurable. □ 

Remark 3.1 We note that the bias and the variance are proportional to e and we emphasize that 
the bias is exactly the generator of diffusion computed for the function F at the starting point Xq, 
the variance is the quadratic variation of diffusion associate to function F at the starting point 
X ; this example show the relation between functional analysis and perturbative theory. 

Finally we can state: 

Theorem 3.1 (impact of uncertainty) The impact of uncertainty on the parameter trans- 
forms a constant into a gaussian distribution of the form 

(3.1) F(a ) - F(a ) + e A[F(a)]\ a=ao + yfe ^T[F(a)}\ a=ao G 

where G is an exogenous independent standard Gaussian variable. 

Remark 3.2 The bias eA\F{a)]\ a=(T0 exists even though the parameter is unbiased; it suffices that 
the function F is non linear. 

This expansion explains the role of the generator and the carre du champ operator (see Bouleau 
[7]); the theoretical image is perturbed due to the uncertainty on the parameter. This effect is 
small, however it produces not only a noise but also it alters the mean. 

In this article we study in particular this shift of mean and we search to reproduce the smile 
effect by means of this shift. 

4 Perturbed Black Scholes Model 

We start with the classical Black Scholes model (see Black et al. [2] and Lamberton et al. [To]): 
let (Q, J- ', P) the historical probability space and B t the associated brownian motion, we suppose 
that the dynamic of the risky asset under historical probability P is given by the following BS 
diffusion without driffi 

dS t = S t a dB t 

In this framework, the price of a European vanilla option is well known (see Lamberton et al. 

[IB]). 

2 We can remark further in this article that the presence of a drift term has an impact otherwise from the 
classical BS model. 
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The B&S model present many advantages, in particular the pricing only depends on volatility 
and we find closed forms for premium and greeks of vanilla options; unluckily the B&S model 
cannot reproduce the market price of call options for all strikes at the same volatility, this effect 
is called smile. 

We propose to consider a perturbation of this model by means of an error structure on volatil- 
ity 

We make three hypotheses: 



1. the real market follows a B&S model with fixed and non perturbed volatility a ; 

2. the trader has to estimate the volatility, so its volatility contains intrinsic inaccurancies, 
we model this ambiguity by means of an error structure; nonetheless we assume that the 
stock price St is not erroneus. We evaluate the impact of the perturbation, generated by 
the trader mishandling, on the profit and loss process used by trader to hedge the vanilla 
option; 

3. the trader knows this perturbation and he wants to modify the option prices to take into 
account the bias induced by the perturbation on volatility. 



4.1 "Mismatch" on trading hedging 

We consider a trader that use an "official" BS asset model in order to hedge vanilla options; he 
uses the market price to determine the "fair" values of parameters in his model (in this case only 
the level of flat volatility <7o) by inversion of pricing formula. 

The trader finds an observed volatility process q, usually known as implied volatility. He 
hedges his portfolio according to his volatility, so the price of an option, that pays a payoff $ is 
Ffa, x, 0). 

We study the profit and loss process associated to the hedging position. 
The profit and loss process at the maturity of a trader that follows the strategy associated 
with his volatility q is given by: 

(4.1) P&L = F(9,, x, 0) + £ S t , t)dS t - $(S T ) 

We make two remarks: 



Remark 4.1 The profit and loss process is stochastic, due to two random sources: 

• First of all, the stochastic "real" model since the trader cannot use the correct hedging 
portfolio. 

• Second, the stochastic process q t , that can depend on a random component independent to 
the brownian motion B t . 

Remark 4.2 The profit and loss process must be studied on historical probability P, therefore the 
presence of a drift on the BS diffusion modifies the second term of equation 4jJ_ Without drift 
this term is a martingale and this fact simplifies the computation. The case of BS model with 
drift will be dealt in a next article. 
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In order to analyze the law of P&L process, it is sufficient to study the expectation on a class 
of regular test functions /i(P&L) and the error on them. 

We will come back on the role and choice of a particular function h in the next subsection. 
We suppose by simplicity that the trader volatility q is a time independent random variable: 

(4.2) q t = a. 

We define an error structure for the volatility a, therefore this volatility admits the following 
expansion: 

ct -> a + eA[a)(a ) + y/ eT[ar](a )Af 

where Af is a standard gaussian variable defined in a space (Q, A, P) independent to Q. Further- 
more, if the volatility has been estimated by means of a statistic on market data, we can specify 
the functional T. In fact, thanks to a result of Bouleau and Chorro [BJ, T is related with the 
inverse of the Fisher information matrix. 

We want to estimate the variance and bias error of E[/i(P&L)]. To perform the calculus, 
we assume that a = <jq is the right value of the random variable in the sense that q = ao and 
P&L(<7o) = 0. We have the following relation for the sharp of volatility: 



Then we can state 



Theorem 4.1 We have the following bias and variance: 



(4.3) 
(4.4) 

where 



TfVo) 
TfVo) 



A[E[h(PLL)]] = h'(0)T^(a ) + -h"(0)T^(a ) 



r [E [h (P&L)]] = [h'(0)] 2 A BS (a 



OF ld 2 F 

— (cr ,x, 0)A[a](a ) + --^(a ,x,0)T[a}(a ] 
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dF 



0"n 
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E 



d 2 F 
dadx 



((T , S t ,t) 



dt\T[a](a ) 



TM(a ) 



and we have the following truncated expansion: 



(4.5) E [h(PSzL)} « eh'(0) Tf 5 (a ) + eh"(0) Tf s (a ) + yj e [h'(0)] 2 A BS Af(0, 1) 
Proof: 

We start with the study of the variance. A computation yields 



(E \h(P&L)]) 



# 



E 



(dF f T d 2 F 

h' (P&L) Ifofa, x, 0) + q^(vo, S s , s)dS, 
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thus the quadratic error is equal to: 



r [E [h (PkL)]] = [h'(0)} 2 lE 



dF 



(a o ,x,0) + 



T 02 



d 2 F 

dadx 



and the second term vanishes since the stock price is a martingale, here the hypothesis on drift 
is crucial. 

The study of the bias is more complicated, we start with the remark that the bias is a linear 
operator: 



A [E [h(PkL)]] = E [A[h{PkL)\] = E 



h\PkL)A [PkL] + ~ti'{PkL)T [PkL] 



We study the two terms separately; for the first, we find the expectation of quadratic error in 
the case q = ctq: 



E [r [PkL]} 



E 



T 



E 



( d 2 F 

a o S s ( ^ZEz( a o, S s , s) ) ds 



V dadx 



T [a] (a ) 



We study the bias operator; we must evaluate the expectation of bias of profit and loss process, 
and we find the following result always in the case q = a : 

dF ld 2 F 
E[A[PkL]] = —(ao, x, 0)A[a]{a ) + - ^j(^o, x, 0)r[a](a ) 

Finally the bias of expectation of a function of profit and loss process: 



OF 1B 2 F 
A[E[h(PkL)]] = hXO){—(a ,x,0)Aia](a ) + - — (a ,x,0)r[a](a ) 



+ 2^(0) 



dF 



+ 



+^0 



T 



E 



S 2 
1 \dodx 



(cr , S t ,t) 



dt\T[a}(a ) 



The proof ends with the truncated expansion that is a consequence of the error theory using 
Dirichlet Forms (see Bouleau [H] and [H]) 

In order to interpret this result in finance, we consider that the trader knows the presence of 
errors in his procedure and wants to neutralize this effect. 

We associate: 



the variance of h(PkL) process to the bid-ask spread of options; 

the bias of h(PkL) process to a shift of prices of options asked by the trader to the buyer. 
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Figure 2: Impact of ambiguity: the Dirac distribution of price X becomes a continuous distribu- 
tion; the mean shifts of e A[X] and the variance is e 



Indeed in the classical theory of financial mathematics we assume that all market securities 
have a single price, with the probability theory language we can associate at any derivative 
securities a Dirac distribution for its price. If we take into account uncertainty on volatility, we 
have found that the price of the contingent claim is not unique but we have many possible prices; 
thus the Dirac distribution changes into a continuous distribution, characterized by a variance 
and a shift of the mean with respect to the previous Dirac distribution (see figure W) ■ 

Theorem 4.2 If the perturbation in the volatility is small, we can neglect orders higher than the 
second, so we always work with Gaussian distributions; the trader must modify his prices in order 
to take into account the two previous effects, namely the variance and the bias, then he fixes a 
supportable risk probability a < 0.5 and accepts to buy the option at the price 

(Bid Premium) = (BS Premium) + e A [E[h(PLL)}} + f [E[h(PLL)}] M a 

where M a is the a-quantile of the reduced normal law. Likewise, the trader accepts to sell the 
option at the price 

(Ask Premium) = (BS Premium) + e A [E[h(PLL)]} + f [E[h(PkL)}} A/i-a 

Remark 4.3 We remark that the two previous prices are symmetric, since Af a + A/"i_ Q = 0; 
therefore the mid-premium is 

(Mid Premium) = (BS Premium) + e A \E[h(PSzL)]] 

We emphasize that with our model we can reproduce a bid-ask spread and we can associate 
its width to the trader's risk aversion (the probability a) and the volatility uncertainty (the term 
T [E[h(PkL)}}). 

In the rest of this article we work directly with the mid premium, but all results represent the 
center of a normal distribution, in order to reproduce the bid and the ask premium we need to 
specify the probability a. 
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We conclude this subsection with a remark. 

Remark 4.4 The presence of the perturbation induces a problem on the completeness of the 
market, the market with perturbed volatility is not complete, since the volatility depends on a second 
random source orthogonal to Q, and the presence of a bid-ask spread is a direct consequence of 
this fact; on the other hand the enforcement that a to be equal 00 cancels the impact of the second 
random source. This apparent contradiction is due to the fact that an argument acts precisely at 
the boundary between complete and incomplete markets. 

4.2 Role and choice of h functions and relative index 

In this subsection we discuss on the choice of function h in equation (14.5p . because function h 
defines the magnitude of correction on prices; this choice becomes simpler since we have to specify 
only the first and second derivative in zero; therefore we can consider that the function h is a 
parabola that passes through the origin. Owing to the two degrees of freedom associated with 
e and a, we can take h'(0) = 1: this is easy to understand from the economical point of view 
because the trader wants to balance his portfolio, i.e. the P&L process. If we look at equation 
(14.51) we find that the choice of h'(0) = 1 defines completely the term of variance. Since the second 
derivative of h has an impact only on the bias and the coefficient T^ s (o"o) is positive, this impact 
is a shift of the mean, as in the following figure. 



h"<0 h"=0 h">0 




Figure 3: Impact of ambiguity: the convexity (resp. concavity) of function h raises (resp. reduces) 
the mean of prices but leaves the variance unchanged. 

We suggest to interpret this impact as an asymmetry of the balance between supply and 
demand. Indeed, if h"(0) = we find that the function h is the identity: this means that the 
trader uses directly the process of profit and loss, and the bias is "neutral" , i.e. that we find the 
same result if we consider the buyer's point of view (it is enough to take minus identity function 
as h). A surplus of the demand of an option with respect to the supply induces a raising of the 
prices: this is the classical case of market where banks sell options and private investors buy. We 
model this perturbation with a positive second derivative for h and we consider that if h"(0) > 
(resp. h"(0) < 0) the demand (resp. supply) exceeds the supply (resp. demand). 

We define the following index of asymmetry of balance between supply and demand: 
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h"(0) 

rs/D - im 

We can identify this index by means of the classical utility theory: if we interpret h as a 
utility function, rs/D is known as the absolute index of h. In the next part we study the bias of 
profit and loss process in the case of a call option; to simplify matters, we suppose that h is the 
identity function i.e. rs/D = 0, but before we must introduce an other index, very important in 
the continuation of this article. 

We concentrate our attention on a problem; after the perturbation of a parameter o"o, we have: 

a ->■ o = a + e A[a] + y/eT[a]J\f with M ~ iV(0, 1) 

However e is generally unknown: the Error Theory via Dirichlet Forms cannot define this 
parameter. In order to deal with this question we propose to renormalize this problem; we 
consider the ratio between the bias and the variance, since the variance is almost surely strictly 
positive, therefore that the dependence on e is cancelled: 

Bias X eA[X] A[X] 



Variance X eT[X] T[X] 
This ratio is not homogeneous, because the generator is linear and the operator "carre du 
champs" is bilinear, so we define a relative index by: 

(4.6) r r (X) = 2^ 

The factor 2 will be justified in section [UJ where we show a relation between Dirichlet forms 
and utility theory, since we can interpret r r (X) as a relative index of an exogenous utility function. 

4.3 Call options case 

We concentrate on call option and we study the bias and its derivatives in order to determine 
some sufficient condition to force the presence of a smile on implied volatility. 

We know the premium of a call option (see Lamberton et al. [IS]) with strike K and spot 
value x, and its hedging strategy: 

C(cr , x, 0) = F(a , x, 0) = xN(d x ) - KM(d 2 ) 
OF 

Delta = —(a ,x,0)=J\f(d 1 ) 
ox 

\nx-\nK+ ^-T / _ 

where d± = -= and d 2 = d\ — aoVT. 

(t VT 

The following results are classical (see 



[<7 , X, U) = XV 1 



da ' y/2% 

d 2 F, , xy/Te-^ 



(4.7) a-2"( a o, x, 0) = — j=did 



d 2 F x e _ ^ d i 

{a , x, 0) 
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Then the bias of the call premium is given by 



(4.8) 



i4 



e 2' 
x—== { A 
V2^ 



d\d 2 



\(T=a 



r 



aVT 



\u=ao 



We can compute the first derivative with respect to the strike: 



dA[C] 



dK 



x e 2 



i ( / 2 

2 1 



Ka VT V27T 
diA[C}\ a=ao 



d x A 



aVT 



di + d 2 - d\d 2 



2aVT 



x e 2 



{d 1 + d 2 )r 



IT 



Ka VT IKolT v^r 

We find that the first derivative vanishes at the forward monej{§ if and only if the bias of call 
vanishes at the same strike. 



8A[C] 



dK 



K=X, <7=0~Q 



e s 
x , <^ A 



2x 



2tt 



A[C}\ 



aVT 



aVT 



K=X, Cf=(7Q 



We can remark that the bias and its first derivative are positive (resp negative) at the money 
if and only if 



(4.9) 



rf s ( aVT 



2anVT 



A 




Ft 


r 


a\ 


FT 



> (resp. < 



We find three cases: 

1. if r r [^J a^T^j < \o^T , then the bias of call and his first derivative are negative at the 



money. 

2. if r r (^J (JqT^ = \(JqT , then the bias of call and his first derivative vanish at the money. 

3. if r r ( (TqTj > \<JqT , then the bias of call price and his first derivative are positive at the 
money. 

Remark 4.5 This bound increases with maturity; if we suppose a constant relative index (CRI) 
we can define a bound on maturity Tuas ■ 

Therefore if we study an option with maturity smaller (resp. greater) than T bias we have that 
the bias associated to the hedging "profit and loss" process and his first derivative are positive 
(resp. negative). 



3 In interest-free case the forward money is for K = x. in this case we have d\ 
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We calculate the second derivative: 



8 2 A[C] 



dK 2 



do 



dA[C] 



KaoVT dK 



x e 2 



if/ 2 

■ 2 a l 



K 2 a 2 T v^F 



.4 



aVT 



| df + 2did 2 -2 r 



and we evaluate the second derivative at the forward money 
8 2 A[C] 



OK 2 



1 e s 













|cr=O"0 



K=x, a=<jQ 

If we force the bias of the call to be convex, we find: 



32 8 



a 2 T 



r 



(T 



(4.10) 



2(7^" 



CT=cr 



A 




/r 


r 




fr 



< 



alT 2 + 4agr + 32 
4<j2T + 16 



6 ((ToV?) . 



J <7=<7o 



Remark 4.6 Previous bound Q(x) is strictly positive, decreasing in [0, Ay/2] and increasing if 
x > AV2, and we find that 9(4^2) = |V2 - 3. 

If the relative index is constant (CRI) we have an always convex bias if r r < 0(4v^2) 

In the previous relations we note that the constraint depends on the volatility by means of 
cumulated volatility ooa/T. For more generality in this study we can assume that the erroneous 
parameter is not the volatility but the cumulated variance f Q a 2 (s)ds that appears in general 
Black & Scholes model when the volatility is deterministic but depends on time. 

Now we study the evolution of slope as a function of maturity at the money 



OA 
df 



<7=<T0 



dA 

dT 



K=x, cr=<xo 



x e 2 
2T7^ 

xe s 



{l + d l d 2 )A 
(4 - olT) A 



rr 



(7 



+ 



Ad\d\ - 3a 2 T - (d 1 + d 2 ) 



o-qVt 



l<r VT 



(a 2 T - 12) T 



-r 



a 



d 2 A 



dKdT 



x 



.4 
e 2 



<T=0-Q 



2a T2 KV2n 



d 2 {l-d\)A 



aVT 



\a=<r 



Ad\d\ - U x o\T + (di + d 2 ){A - 9d t d 2 - d\ 
+ d\ ^= 1 



8a VT 



aVT 



I a=(TQ 



d 2 A 



dKdT 



e 8 



K=X, (T=(7Q 



16TV2T 



(a 2 T -A) A 



aVT 



\<j=a 



aoVT 



alT-Yl. 



a\ff 



\(T=a 



to find a slope that increases with increasing maturity we have to impose: 
(4.11) 



r r ( aVT) \ a=ao > 



alT alT - 12 , ~ 

-9— -r o — with olT < A 

A A — a 2 T 
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We study the evolution of smile as a function of maturity. 



3 A 



dK 2 dT 



"0 T 

e s 



16 



xa 2 T 2 y/2^ 



rr 4 T 2 



K=X, 0~=(TQ 

This term is positive if and only if: 



32 



aVT 



a 2 T(a 2 T- 4) 2 + 128 r 



(7> 



\cr=ao 



256 



aVT 



a=cr 



(4.12) 



4.4 Example 



r r [aVT] \ a=ao > 



1 o%T (alT - 4) 2 + I2i 
4 



1 + a^T 2 



Now we study a particular case, in order to reproduce the smile of volatility present in market 
data. 

Theorem 4.3 We fix the relative index at 



(4.13) 



r? S (aVT)\ c 



a 2 T 



This choice fix the values of the bias and its derivatives at the money, in particular this choice 
force the bias and its first derivative to be zero at the money (see equation \4-9\ ); the second 
derivative becomes positive at the money thanks to equation \4-10\ for any maturity T. Therefore 
the bias is strictly convex around the money and vanishes at the money, therefore it is positive in 
a neighbourhood of the money. Now if the bias vanishes at the money the implied ATM volatility 
is o"o, but, since the bias is positive around, the implied volatility becomes greater than o"o around 
the money^; Finally we have reproduced a smile effect around the money. 



4.5 Dupire formula and Implicit Local Volatility Model 



In this section we want to specify the Local Volatility Model equivalent to Perturbed Black Scholes 
Model. We know that the knowledge of prices of options for all strikes and maturities defines a 
single local volatility model that reproduces these prices; the Dupire formula (see Dupire [TO] ) 
defines the local volatility function: 



(4.14) 



<p(T, K) 



ac 

dT 



1 K1 d 2 C 
2 X1 dK 2 



Our model is a perturbation of BS model, so we can consider the following expansion: 

BS Term Perturbation 



C(cr, x, K, t, T) 



dC 

— (a, x, K, t, T) = t^(^o, x, K, t, T) 



C(a , x, K, t, T) + 



dC 



eA[C}(a, x, K, t, T)\ 
dA[C] 



dT 



d 2 C 
dK 2 



a, x, K, t, T) 



d 2 C 
8K 2 



(a , x, K, t, T)+ e 



dT 

8 2 A[C] 
dK 2 



(a, x, K, t, T) 
a", x, K, t, T) 



a=a 



4 We recall that the vega is positive for call options. 
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Figure 4: Representation of the smile in PBS model. 



But in fact, if e vanishes, the model is a Black Scholes model with volatility a®. Thanks to 
equation (14.71) . we can rewrite: 



o? m JT, K) « al 



1 K2 d 2 C 
2 ly dK 2 



dA[C] 



dT 



2_2 



d 2 A[C] 



a=(T 



dK 2 



(7=170 



Theorem 4.4 The local volatility model equivalent to the PBS model for vanilla options, has the 
following local volatility function: 



(4-15) o-L P (T, K)^al{l + - 



A 



a VT 



a 2 T + 2 



4m(f) : 



alT 



Remark 4.7 Local volatility o~(T, K) has a minimum at forward money K = x, and present a 
logarithmic behavior as K approaches zero and infinity. 

We must preserve the positivity of the square of volatility, so we fix the following constraint: 



alT + 2-2r r (aVf) L =CT() < <^ e 



r 



a 



T 
ofT 



5 Extension and relation with literature 

We have proved that the PBS can reproduce at the same time the bid-ask spread and the volatility 
smile; but in literature many authors (see Perignon et al. [T7j, Renault et al. [IB] and Rubinstein 
[T§] ) have remarked that, generally, the volatility presents a skewed structure (the graph of implied 
volatility is downward sloping); besides in this article we have limited the study of the PBS model 
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Figure 5: Examples of local BS and PBS variances (square of volatilities) for ATM = 100, a\^T = 
30%, e = 0.01, T[aVT] = a 2 T and r r = 0. 



at the martingale case, when the extra returns of the stock are zero; however, a simple argument 
of risk aversion induce leads us to suppose that the parameter \x (the extra-returns term in the 
BS model) must be positive. In a next paper we will show that the extra returns term has an 
impact in PBS model, contrary to BS model, and we can use this effect to generate a slope in 
implied volatility. We want emphasize that the PBS model uses a perturbative approach, i.e. we 
start with a simple model (the Black Sholes model) and we adjust the model at the market data 
through a perturbation of the principal parameter (the volatility). 

Clearly in literature some authors introduce a perturbative approach in finance, we recall the 
papers of Hagan et al. [H] and the book of Fouque et al. [IT] ; our approach is however different, 
it is a probabilistic approach and it can relate the bid-ask spread and the volatility smile through 
the simple economic argument of the existence of uncertainty in the market. 

The study of SABR model, a stochastic volatility model introduced by Hagan et al., is based 
on a small volatility of volatility expansion. Now the principal difference between our perturbation 
method and the Hagan's one is that our model is based on a probabilistic point of view, therefore 
our perturbative approach takes into account the second order derivative and, in particular, it can 
estimate the bias induced by non-linearity. In the paper of Hagan et al. the approach is based on 
analysis, so that they implicitly suppose that an exact value for the volatility and, by consequence, 
for the option price exists, their perturbative approach is therefore a first order formulation; they 
cannot then justify the presence of a bid-ask spread directly through their model. 

Moreover Hagan et al. start with a more complex model (a stochastic volatility model with a 
volatility driven by an brownian diffusion with an independent component) and they simplify their 
model with the perturbative approach in order to find a closed form, therefore in SABR model 
the perturbative approach is a computation tool. Our approach allows to justify the perturbative 
approach with an economics reason, the problem for the traders to estimate the value of volatility. 
Therefore in our method the perturbation is naturally a consequence of uncertainty on volatility, 
and it has a "physics" nature; finally our model start by the classic, well-known and accepted 
Black Scholes model, in one phrase the SABR model is a downward perturbative approach the 
PBS is an upward one. The principal advantage of the SABR model with respect to PBS model 
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is that the SABR model can predict the dynamics of the smile when the forward price of stock 
changes; this capacity is related with the parameter f3, and has no relation with the perturbative 
approach. 

In the book of Fouque, Papanicolau and Sircar the authors study a model based on the Black- 
Scholes model with a stochastic volatility with a fast mean reversion, in this case the perturbative 
approach is based on the presence of two time scales the mean reversion time and the maturity of 
options, the first is very small compared the other time scaled, therefore, in their book, Fouque 
et al. introduce an expansion of the solution of their prices PDE (based on a hierarchical PDE 
system). Therefore, in their book, Fouque et al. use a perturbative approach with a "physical" 
nature, i.e. the presence of a two time scales; but their approach is based on analysis and it 
is not (directly) related with an uncertainty on volatility; finally the Fouque's approach cannot 
reproduce a bid-ask spread. 



6 Risk Aversion 

In this section we make some recalls on the theory of utility functions and we show a connection 
with Error Theory using Dirichlet Forms. We consider a utility function U(x) : 1R — > K and 
U{x) G C 2 : let p be defined by the following relation 



(6.1) E [U(X)\ = U (E [X] - p) . 
We find, in the case of a small variance, that: 

(6.2) ^■■( EW W 

v ; H 2 U'{E[X)) 2 aK 1 

In a similar way we define p so that: 

E[U(X)} = U {E[X}(1 - p)} : and we find 

We can name the three objects: 

1. p is the risk price; 

2. r a (X) is the absolute index of aversion at the wealth X; 

3. r r (X) is the relative index of aversion at the wealth X. 

We write the relations of error and bias, given an error structure on X: 

A[U{X)\ = U'{X)A[X) + ^U"{X)r[X] 
T[U(X)] = (U'(X)fT[X] 

We observe that the bias of U(X), A[U(X)] is zero if and only if A[X] = p; in this case we 
find 
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A\x] = ±r[X\ 

(64) ffl = 

We have found a relation between the utility function theory and the error calculus using 
Dirichlet Forms. We suppose that all traders buy and sell according to their risk aversion, and 
this aversion is represented via a utility function U(x); then the utility function is defined by its 
relative index of aversion. 

We make an hypothesis: 

Hypothesis (*): for a trader the bias of the utility of a traded wealth vanishes. 
We can interpret this hypothesis from an economics point of view in two ways: 

1. the utility function of trader, supposed to be known, is like a lens traders look the market 
through. Traders don't add any effect to balance their aversion; 

2. the vector (X, A, T) is supposed to be known, we can define the utility function of a trader 
as function U (x) that cancels the bias of U(X) where X is the considered wealth. 

Under hypothesis (*) we have two relations between the utility function and Dirichlet Forms: 

1. A[X] = p(X), where p is the risk price; 

2. r r (X) = 2^ 

Remark 6.1 Thanks to relation \4-l'J\ we can define the class of utility function that preserve 



vanishing bias and slope at the money. 



X 



U"{X) X 2 

U'(x) T 

U'(X) = 
U(X) = 

Where Af(X) is the distribution function of the normal law. This utility function is concave 
if the wealth is positive and convex otherwise. 



7 Conclusion and Economics Interpretation 

In this paper we have studied the impact of a perturbation on volatility in Black-Scholes model in 
absence of drift and term structure; in particular we have dealt with the problem of call hedging. 

We have proposed a new model for option pricing, called Perturbed Black-Scholes model; the 
basic idea is take into account the effect of uncertainty of volatility value in order to reproduce, 
at the same time, the spread bid-ask and the smile on implied volatility; the mainly advantage 
of this model is that it is based on the classic Black Scholes model and the price of an option in 
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PBS model is the BS price plus a small perturbation that depends only on the greeks founded 
with BS formula, therefore the computation of PBS price is given by a closed form. 

The PBS model depends on four parameters, naturally on the volatility of stock; but also on 
the variance of the estimated volatility er[cr](a"o), on a relative index r r (a), that represents the 
ratio between the bias and the variance of the estimated volatility, and, finally, contrary to Black 
Scholes model, on the drift [i that represents the extra returns of stock; the impact of the drift 
rate will be studies in a next paper. 

In particular, if we set the relative index to be equal to T , we have proved that the implied 
volatility present a smile around the money. 

Finally we have defined a Local Volatility Model equivalent to Perturbed Black Scholes Model 
as far as vanilla options are concerned; the related local volatility function is defined by Dupire 
formula: 



a\T, K) naUl 



.4 



a 2 T 



4m(f) : 



aVT 
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